A POSSIBLE COUNTEREXAMPLE TO WELLPOSEDNESS OF 
ENTROPY SOLUTIONS AND TO GODUNOV SCHEME 
CONVERGENCE 

VOLKER ELLING 



Abstract. A particular case of initial data for the two-dimensional Euler 
equations is studied numerically. The results show that the Godunov method 
does not always converge to the physical solution, at least not on feasible grids. 
Moreover, they suggest that entropy solutions (in the weak entropy inequality 
sense) are not well-posed. 



1. Introduction 

Consider the Cauchy problem for a system of hyperbolic conservation laws, 

(1) ^ + V- (/(«))= 0, 

(2) u(0,-)=u , 

where u = u(t, x) : R+ +1 := (0, oojxl^Fc R m is the desired solution (P 
the set of physically reasonable values), / = (/*), /* : P — > M. m , the (smooth) flux 
function, uq : M. d — > P initial data. Here and in the sequel "V, A, ■" are meant with 
respect to x. 

An important example of hyperbolic systems of conservation laws are the (non- 
isentropic) compressible Euler equations: 

Pt + V • ( P v) = 0, 
{pv l ) t + V • { P v l v) +p Xl =0 (i = 1, . . . , d), 

(3) (pe) t + V-((pe+p)v) = 0. 

Here, p is density, v = (v l ) velocity, e specific energy, which decomposes into 
(4) 

the first summand is specific kinetic energy, q is specific internal energy. The 
pressure is a function of p, q; a common choice is the polytropic pressure law 



(5) p = (7 - l)pq 
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(1 < 7 < |; for air, 7 = |). The set of admissible values is 

P={q>0, p>0}. 

It is well-known that |QJ and J2J need not have a global smooth solution, even 
if the initial data uq is smooth. For this reason, one has to study weak solutions, 
defined as functions u £ i 1 1 oc (R ( | +1 ; P) that satisfy 



(6) II u— + f(u) ■ V</> dx dt = I u a (x)(j)(0,x) dx, 



for all test functions (f> £ C^°(IR^ +1 ). Moreover, there can be more than one weak 
solution, so it is necessary to impose an additional condition, called entropy condi- 
tion, to single out a unique weak solution (the entropy solution). 

One definition of entropy solutions is the vanishing viscosity (VV) definition; it 
requires that u is the limit of the sequence (u £ ) e >o of solutions of 

(7) ^! + v.(/V)) = eA« e mRf 1 , 

(8) u e (0,-) = u on{0}xM d . 

The limit is taken in some suitable topology, usually as a boundedly almost every- 
where limit. We call such a function u a VV solution. 

Another definition uses entropy / entropy flux (EEF) pairs (rj, tp), where rj : P — > K 
is a smooth strictly convex function, called entropy, whereas "0 = (ip 1 , ■ ■ ■ , ip d )' with 
smooth : P —> M is called entropy flux; r\ and ip are required to satisfy 

( 9 ) d^ = 22^^ (i = l,..;d,a = l,...,m). 

(3=1 

By multiplying from the left with r]'(u e ) and using ©, one obtains 



(10) ^ + _ eA(n o u<) _ eE „, V) g - s eA( , 



(here, we used that 7/ is convex). Upon multiplying the last equation with a non- 
negative test function </> and integrating by parts, this yields 

(11) 



r](u € )-^- + ip(u e ) -V4>dxdt<el / 7](u e )A0 dx dt + / 7?(u ) ckr. 



o 



If, as assumed above, (u c ) — > it boundedly almost everywhere, then implies 



(12) - / / Tl{u)jf- + i'(u) dx dt< / 77(1*0) dx. 

Functions m that satisfy (|12|) for aZZ EEF flux pairs are called EEF solutions (of 
Q). As we have shown, VV solutions are necessarily EEF solutions. 

In the literature, the term entropy solution is used to refer either to EEF or to 
VV solutions, often without explicit mention, because it has been assumed that the 
two definitions arc equivalent for the Euler equations and many other physically 
relevant systems (see |Ser99| p. 101, |Daf00| p. 49, |GR96j p. 32; sec the discussion 
in Section for verified special cases). However, for the purposes of this paper it 
is necessary to distinguish the two notions, as we will discuss a possible numerical 
counterexample to their equivalence. 
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The (gas-dynamic) specific entropy s is defined as 

(13) s = logq+ (1 -7)log/o; 

(14) rj := -ps, if,' 1 := -psv 1 

provides an EEF pair for the Euler equations. 

A common simplification is to assume that s is constant in space and time. This 
yields the isentropic Euler equations 

p t + V • (pv) = 0, 

(15) {pv% + V-(pv%+p Xi =0 (i = l,...,d) 
with 

(16) P (p) = pi. 

In this case, P = {p > 0}. An EEF pair is provided by the specific energy e, 

\v\ 2 p<~ x 



2 7-1 



with 



r) := pe, V* : = (P e +p)v l - 

It is cumbersome to verify the EEF condition 1)12(1 directly, not to mention the 
VV condition. There are easier criteria for piecewise smooth functions, which we 
define in the following customized way: 

Definition 1. (1) A point (t, x) E is called point of smoothness if u is 

C°° in a small neighbourhood of (t, x). 

(2) A point (t, x) S R + +1 is called point of piecewise smoothness of u if there is a 
C°° diffcomorphism $ of a ball V around in R d+1 onto a neighbourhood B 
of (t, x) = <f>(0) so that u off) is C°° on 5_ and on B + (where B± := 3>(V±), 
V± := {y eV :yi^ 0}); for later use, let S be the surface $(^n({0} xR d )), 
n = (n*, n) G a unit normal to S in (i, x) pointing into £?+; let u+, u_ 
be the one-sided limits of u in (i, x) within 5_ resp. -B+). We also require 
n ^ 0. 

(3) u is called piecewise smooth if there is a set A of d-dimcnsional Hausdorff 
measure so that all points in — N are points of piecewise smoothness. 

Proposition 1. Let u be piecewise smooth, u is an EEF solution of Q if and 
only if 

(1) it is a (classical) solution of Q in each point of smoothness, 

(2) u(t, ■) — * wo in L\ oc as t [ 0, and 

(3) in each point (t, x) of piecewise smoothness it satisfies the Rankine-Hugoniot 
conditions 

(17) (u+ - + (/(«+) - /(«_)) • n = 
and (7or aH SEF pairs (?y, j 

(18) (r)(u+) - r?(w-))n* + ($u+) - $u_)) • n < 0. 

Proposition is well-known (see, for example, Section 11.1.1 in |Eva98| ). as is 
the following property: 
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Proposition 2. For the Euler equations lO resp. ()15|l (with polytropic gas law © 
resp. Ijl6(l ), Ijl8|l is equivalent to the simpler condition that the normal velocity does 
not increase across discontinuities: 

(v+ — V-) ■ n < 0. 

The Cauchy problem for the Euler equations has several important symmetry 
properties, including the following: 

Proposition3. Letu — (p,v' ,q)' be a weak solution for initial data uq — (po, v' , qo)' ■ 

(1) Change of inertial frame: For allw G M. d , (p(x+wt),(v(x+wt,t)—uiy,q(x+ 
wt))' is a weak solution for the same initial data uq. 

(2) Self- similarity: a function f : — > M m is called self-similar if f(rt, rx) = 
f(t, x) for all r > 0; same for functions on M d . // the initial data is self- 
similar, then for any r > 0, u(rx, rt) is a weak solution for the same initial 
data uq. 

These symmetries remain true after replacing "weak" by "VV" or "EEF" . Analo- 
gous symmetries hold for the isentropic case. 

2. Example and numerical results 

Consider the following set uq of initial data for © with d = 2 (see Figure ^1 : 
the data is symmetric under reflection across the x-axis and constant in each of 
four cones centered in the origin (in particular, constant along rays starting in the 
origin). In the origin, two shocks emanate into the first and fourth quadrant; the 
area on the left is supersonic inflow (parallel to the x-axis); the two areas on the 
other side of the shocks are denser and hotter gas, moving parallel to the contact 
discontinuities (see [CF48j Chapter IV C on choosing pre- and post-shock values 
that satisfy the Rankine-Hugoniot conditions; we choose the ones that yield the 
weaker shock). The gas in the stagnation area (enclosed by the contact disconti- 
nuities) has the same pressure as the post-shock gas on the other side, but velocity 
v = 0. It is easy to check, using Propositions ^ an d El that the steady solution 
u(t,x) = uo(x) is an EEF solution of J2Jl resp. (|I5J) . Henceforth we refer to it as 
Solution T (for theoretical). 

However, instead of Solution T, numerical calculations produce the markedly 
different result in Figurc|2]which we call Solution N (for numerical) in the sequel (of 
course it is not known to be an exact solution). The numerical domain in Figure|21is 
indicated as the dotted quadrilateral in Figure^ Figurc[21was computed as follows: 
adaptive refinement was used to achieve better resolution at same computational 
cost. To reduce numerical viscosity the grid was chosen so that near the right 
domain boundary the edges are aligned with the contact discontinuity and the 
shock. In order to capture self-similarity, the computations were done for a grid 
with moving vertices with coordinates x = i£ (£ has the dimension of a velocity; 
its components are called similarity coordinates). The moving-edge modifications 
discussed in E1100| Section 2.1.6] and (Eu05l Chapter 4] were used (the essential 
idea is to compute numerical fluxes across a moving edge by transforming to a 
steady edge, using invariance under change of inertial frame (Proposition |3Jl, and 
to apply an arbitrary approximate Riemann solver to the transformed problem). 
The domain boundaries were chosen so that small perturbations on them propagate 
into the domain (£ • n (n outer unit normal) in each boundary point is larger than 
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Figure 1. Solution T: this initial data is also a steady and self- 
similar solution for the compressible Euler equations in 2D 

the maximum of \v\ + c in the domain). This allows to prescribe all components of 
the fluxes on the boundary. 

Experiments with various modifications were made: changing the numerical 
scheme (the experiments were repeated for the Godunov scheme |God59| , the Oshcr- 
Solomon scheme |OS82| . the ENO-RF scheme |SQ89| . and a second-order MUSCL 
code based on the first-order ENO-RF scheme), adding more numerical dissipation, 
refining uniformly rather than adaptively, using a Cartesian grid including origin 
and lower half-plane, or calculating in space rather than similarity coordinates. 
None of these modifications change the numerical results significantly; in all cases, 
the numerical results converge to the same Solution N. 

Solution N appears to be self-similar (i.e. steady in similarity coordinates), but 
it is strongly unsteady, so it is clearly different from Solution T. 

3. Conclusions about numerical methods 

While the discrepancy between Solutions N and T opens many new problems, 
we can already draw one definite conclusion. 

If we assume that Solution T is the correct solution, many 1 popular numerical 
schemes fail to converge to physical solutions. Although it cannot be ruled out that 
they ultimately converge to Solution T as the numerical grid becomes infinitely 
fine, they approach Solution N for computationally accessible grids — which is all 
that matters for practical purposes. 

On the other hand, if we assume that Solution N is the correct solution, there is a 
trivial theoretical example of misconvergence: consider the (scmidiscrete) Godunov 
scheme on grids whose edges are exactly aligned with the discontinuities of Solution 



in the sense of: every scheme that was tested 
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FIGURE 2. Solution N: (Rotate clockwise by 90° to align with 
dotted area in Figure Q]) Each square in the coordinate grid cor- 
responds to a 100 m/s x 100 m/s square in the £ plane. Origin 
marked by diamond (lower left corner). Plotted: horizontal veloc- 
ity. Godunov scheme for iscntropic Euler equations; data: 7 = 1.4, 
a = 10°; inflow: p = 1.19 kg/m 3 , v = 1000 m/s, T = 20°C. The 
solution differs significantly from Figure ^ Results for nonisen- 
tropic Euler equations or other numerical schemes are similar. 
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T (see Figure P): in exact arithmetic it would have Solution T as steady state on 
every grid. 

In either case — even if Solution N is correct, which would be less catastrophic 
for numerical analysis wc have to conclude that discrete entropy inequalities are 
not sufficient to avoid convergence to unphysical solutions on feasible grids. Hence 
they lose a bit of their value as design principles for numerical schemes, although 
they are still useful as easy-to-check necessary conditions that arc sufficient for 
scalar conservation laws and (probably) ID systems (as supported by the recent 
work on small total variation solutions described in Section 

Although many reports of deficiencies of various numerical schemes have been 
published, the clear case of failure observed here has no precedent. 

4. Theoretical interpretation 

It remains to discuss which of Solution T and Solution N is the physical one 
and, if Solution T is correct, what causes numerical schemes to produce Solution 
N. There are three possible explanations (which are not mutually exclusive): 

(1) either Solution N is an example of failure of numerical methods, or 

(2) EEF solutions are not stable (in the sense of continuous dependence on 
initial data), or 

(3) EEF solutions are not unique. 

4.1. Breakdown of numerical methods. It has already been shown in Section 
13 that the Godunov scheme is flawed, in the sense that it can fail to converge to 
the physical solution on feasible grids. Hence it is natural to suspect that Solution 
N is a numerical artifact that does not correspond to a seccond EEF solution (or 
any other type of solution of the Euler equations). 

Solution T is steady and self-similar. "Steadyness" is a non-generic property 
that is usually not inherited by finite-accuracy numerical solutions (for example 
for a Ricmann problem that is solved exactly by a single shock, most numerical 
schemes produce small additional waves and a slightly different shock). However, 
one would expect numerical approximations to be at least almost steady, unlike 
Solution N. 

In a single space dimension, the conservation property of numerical schemes often 
guarantees accurate shock locations, even if the overall accuracy of the scheme is 
poor. On the other hand, in two or more dimensions numerical imprecision can 
significantly change the shape and location of shocks. This may be the cause 
of Solution N. For example, the upwards deflection of the incoming flow by the 
high-pressure area in front of the stagnation region could be weaker in numerical 
calculations than in Solution T; the additional pressure would cause the stagnation 
region to collapse. 

However, in this case the numerical results would depend strongly on the choice 
of numerical method, mesh width and other parameters. This is not observed; 
rather, all choices produce essentially the same results. 

4.2. Instability. A second explanation is instability (in the sense of lack of con- 
tinuous dependence on the initial data). It is possible that Solution N is an ap- 
proximation to an unsteady EEF solution T' that results from a slight perturbation 
of solution T at initial time (such perturbations are inevitable in most numerical 
computations due to inexact arithmetic, discretization error, artificial viscosity etc.) 
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Since Solution N is produced (up to minor differences) for any "perturbation" (i.e. 
for any choice of mesh, numerical method and parameters), it would indicate that 
Solution T constitutes a set of initial values for which the Euler equations are not 
stable. 

On inspection in similarity coordinates, it appears that the numerical solutions 
are bounded, converge quickly to Solution N and remain steady (many orders of 
magnitude of time have been observed), so Solution T" would have to be at least 
approximately self-similar. If we assume it is asymptotically self-similar for large 
time, then the following theorem implies that the asymptote is an EEF solution: 

Theorem 1. Let u € L rx '{M. d ^ 1 ) be an EEF solution of @. A ssume that u is 
asymptotically self-similar (see Definition^), then its asymptotic limit w (a self- 
similar function) is an EEF solution as well. 

(The proof of Theorem ^ and an analogous result for steady solutions are pre- 
sented in the Appendix.) The asymptote would assume the same initial data as 
Solution T, but would have to be different from it (by closeness to Solution N). 
This would already imply the third explanation (nonuniqucness of EEF solution) . 

To avoid that, it is necessary to assume either that Solution T' is approximately, 
but not asymptotically self-similar — for example it might oscillate periodically 
around some self-similar function without approaching it — , or that the self-similar 
asymptote has data at infinity that does not match the initial data (Solution T). 
Either of these cases would be revealed by a sufficiently fine numerical grid: the 
numerical computation on that grid would refuse to converge to a steady state for 
the given boundary data. But the mesh used to compute Solution N is already 
rather fine, as can be seen from the curved shocks in Figure there is no obvious 
reason why an even better grid is required. 

4.3. Nonuniqueness of EEF solutions. The third explanation is nonuniqucness: 
Solution N corresponds to an EEF solution that assumes the same initial data as 
Solution T. 

In this context, the following peculiarity is important: according to Proposition 
^ the EEF condition 1)12(1 is "insensitive" to sets with (d— l)-dimensional Hausdorff 
measure (such as a single point, for d = 2); e.g. if (|12J) is satisfied for <fi £ 

C~(R^ +1 - {0}), it is satisfied for all $ E C~(R£ +1 ). In verifying that Solution 
T is an EEF solution, we may ignore the rather singular wedge tip in the origin 
of Figure ^ It is counterintuitive that conditions for the physical correctness of 
solutions may ignore such singularities. 

It seems unlikely that the solution in Figurc^is stable under small perturbations 
at the origin (such as perturbations from viscous terms in the VV limit) . Note that 
a planar shock, with inflow state on one side and stagnation area state on the 
other side, would not be steady but move into the stagnation zone quickly — it 
seems unlikely that the example data, which has less mass and energy and more 
x-momentum in the {x > 0} halfplane, would yield a steady pattern in the origin 
(however, the "maximum principle" implicit in this argument is merely heuristic 
and may be wrong in some instances). 

But it is the following observation that provides the strongest argument for 
nonuniqucness: the Lax-W endroff theorem (see |LW60j : see also |GR96| . KRW96 
and most generally |E1103) for Lax-Wendroff-type theorems for irregular grids) states 
that if a numerical scheme is consistent and satisfies a discrete entropy inequality 
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(see [HHL76I IMQ791 IOC84I ITadMl ITad87l IOT88j ) . then the limit of a boundedly 
almost everywhere converging sequence of numerical solutions is an EEF solution. 
(Note that we prescribe the full numerical flux on the boundaries, so the boundary 
conditions are analogous to an initial condition and can be treated with a straight- 
forward modification of the Lax-Wcndroff theorem for initial- value problems.) The 
Godunov scheme, used to compute Figure [3 is the standard example for a consis- 
tent scheme that satisfies all discrete entropy inequalities. Our numerical solutions 
do, on inspection, appear to converge quickly; this would imply that Solution N 
corresponds to an EEF solution. 



5. Related work 

For multidimensional scalar (m = 1) conservation laws with arbitrary /, |Kru70| 
(generalizing earlier work) shows that a global EEF solution exists, is unique, sat- 
isfies the VV condition as well, and is stable under L 1 perturbations of the initial 
data. 

Gli6£] provides a famous existence proof for strictly hyperbolic systems with 
genuinely nonlinear fields and initial data with small total variation; the interaction 
functionals constructed in this paper are a crucial ingredient for all subsequent 
work. |Liu81| extends the result to systems with linearly and some nonlinearly 
degenerate fields. |BCP00| constructed the Standard Riemann Semigroup (SRS), 
an L 1 -stable semigroup of EEF solutions for initial data with small total variation, 
for strictly hyperbolic systems with genuinely nonlinear or linearly degenerate fields 
(see also |LY99| L |BL97j showed that EEF solutions to ID systems are unique and 
coincide with the SRS solutions, under certain smoothness assumptions including 
small total variation (see also |BG99| ). |BB01j prove that for small TV initial data 
and strictly hyperbolic (but otherwise arbitrary) systems VV solutions exist and are 
stable under L 1 perturbations of the initial data, so for the class of solutions that 
are subject both to |BL97| and to |BB01) . EEF and VV solutions are equivalent. 

On the other hand, an EEF pair (77, ijj) has to satisfy the condition © which 
is an overdetermined problem for m > 3, so for some systems no EEF pairs exist 
and the EEF condition is void. However, EEF pairs do exist for most physically 
relevant systems, even those with m > 3. More seriously, for certain 2x2 systems 
(with nonlinear degenerate fields) |OL81| construct a single weak shock that is an 
EEF solution but does not satisfy the Liu entropy condition (see |Liu74l lLiu75j ). 
By |BB01j , there must be a VV solution (for the same initial data) that satisfies the 
Liu entropy condition as well — so it cannot be the aforementioned weak shock. 
Therefore the example in CL81 also constitutes an example of a nonunique EEF 
solution, albeit for an "artificial" system with nonlinear degeneracy. 

|Hop67| proposes the EEF condition for scalar conservation laws (m = 1), proves 
that it is implied by the VV condition under some circumstances and notes that 
there is a large set of convex entropies. Apparently independently, |Kru70j obtained 
analogous results for systems. |Lax71| contains the first use of the term "entropy 
condition" for the EEF condition. Various forms of the EEF condition had been 
known and in use for special systems such as the Euler equations for a long time 
(e.g. by the name of Clausius-Duhem inequality), especially as shock relations; how- 
ever, the above references seem to be the first to define the general notion of strictly 
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convex EEF pairs, to propose the EEF condition as a mathematical tool for arbi- 
trary systems of conservation laws and to formulate it in the weak form Ijl2(l rather 
than the special case (I18|) . 

|LZY98j provide an analytical and numerical discussion of 2D Riemann problems 
for various systems including the Euler equations. However, they focus on data 
constant in each of the four quadrants, so Solution T is not covered. 

6. Conjectures and final remarks 

The results demonstrate that 

the Godunov method does not always converge to the physical solution on feasible 
grids. 

Moreover, they suggest the following conjecture: 

EEF solutions to the multidimensional Euler equations are not always unique. 
If this conjecture is true, it would have far-reaching consequences. The EEF con- 
dition would not be sufficient as a selection principle for physical/unique solutions, 
except in special cases like the ones described in Section |SJ It would be neces- 
sary to find ways to use the cumbersome VV condition or to discover new entropy 
conditions. 

Although the numerical results support the conjecture unambiguously, the ques- 
tion is so important that a rigorous proof is highly desirable. However, since the 
initial data has large vorticity at the contact discontinuity, it seems difficult to con- 
struct (or to prove results about) exact solutions. One possible line of attack is to 
derive novel entropy conditions from the VV condition and to check whether they 
are violated by the steady solution in Figure ^ 

In any case, this paper motivates the investigation of multidimensional Riemann 
problems for systems; these appear to be very difficult and exhibit a large variety of 
phenomena (see |LL98llL%Y98] ). This goal requires techniques for proving existence 
of smooth steady or self-similar solutions to boundary-value problems for systems 
of nonlinear hyperbolic conservation laws; while there arc classical methods for 
smooth solutions in hyperbolic regions, work on tools for the elliptic and mixed 
case has begun only recently (see e.g. [EL05]). 

Appendix: asymptotically steady and self-similar weak solutions 

Remark: in Theorem^ and in the following statements, 

f(t, -)^g in Ll c (n) 

as t I resp. t f oo is to be understood as: for all e > and K <s Q, there is a 
T = T(e) > so that for almost all < t < T resp. t > T, 

\\f(t,-)-g\\ LHK) <e. 

Lemma 1. Let u G L°°(R+ +1 ), u Q E L°°{R d ). If 

(1) u(t, •) -> u Q in Lj oc (R d ) as 1 10, and 

(2) u satisfies ® for all <f) S C c °°((0, oo) x M. d ), 

then u satisfies © for all 4> e C c °°([0, oo) x R d ). 

Proof. Let 9 e C°°[0, oo) so that 9 = 1 on [0, 1], 9 = on [2, oo). For any T > 0, 
define 9 T (t) := 9{T- 1 t). Note that \9 T \ = O(l), \9' T \ = C^T" 1 ) (as T | 0). For 



WELLPOSEDNESS OF ENTROPY SOLUTIONS 



11 



;iven <j> 6 C c °°([0,oc) x R d ), split (f>i(t,x) := 9 T {t)(/)(t,x) and <f> 2 = (f> - fa. 
ufa + f(u) ■ V0 d.T 



2T 

W>it + /(it) ■ V0i dx dt + j I ufat + f(u) ■ V02 dx dt. 

Since fa £ C C ([T, oo) x K d ), the second summand vanishes by assumption. The 
first summand equals 

= o(t sup ||w(t, ■) - w ||i • + 1)) + / [ u fa t dxdt 

\ 0<t<2T / Jo JR d 

+ 0(T(||«o||i+ sup ||«(t,.)-uo||i 

\ 0<t<2T 



0[ sup \\u(t, ■) — t*o||i ) — / u o (x)(j}(0,x) dt 

\0<t<2T 



+ T(|KHi+ sup ||«(t,.)-uo||i 

\ 0<t<2T 

On taking T J, 0, all O terms vanish; hence u satisfies (0. □ 

Remark: the converse of Lemma ^ (which is not needed) is not immediate be- 
cause u(t, •) — u may be large for some t as long as the set of such t has small 
measure near 0. 

Definition 2. (1) A function u £ L 1 1 oc (R^ +1 ; IR m ) is called asymptotically self- 
similar if there is a function w : WL d -> K m so that 

u{t,t~ l -)^w in L\ oc (R d ). 

(2) u is called self-similar if, for some w, u(t,t ■) = w for almost all t > 0. 

Proof, (of Theorem^! By Lemma^ to show that w is a weak solution it is sufficient 
to check that 



(19) r I " «, (f ) ^(t, + /(«, (f )) . V^(t, : 



dx dt = 

for all 4> £ C£°(]0, oo)). The essential idea is to scale coordinates to shift the support 
of (f> into a large-i region and to use asymptotic convergence. 

Let < t\ < t2 be such that supp</> C [^1,^2] x K d - Let e > be arbitrary, set 
T = T(e) as in Definition [3 The change of coordinates t = ^r, x = ^£ changes 
the left-hand side of l|19|) into 

f) St H I l f ) + /(M(T ' 0) ' ^ H * " T 



(20) 



■T -eT d 
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where O is with respect to e — > oo. Note that the support of the scaled <p is m 
[T, oo) x M. d . Also, the assumption that u is bounded is essential here. The second 
summand on the right-hand side equals 

JM d tl \ 1 1 J T 

Since u is assumed to be a weak solution, this term vanishes. Taking e J, in i|2U|) 
yields ifl^l . 

For the proof of the EEF part, replace u,w by rj(u), rj(w) and f(u),f(w) by 
il){u),ip{w) above. □ 

The same results as for self-similar weak solutions can be obtained for steady 
solutions: 

Definition 3. (1) u e Lj oc (l| +1 ;I ra ) is called steady if, for some w : R d -> 
R m , u(t, ■) = w for almost all t > 0. 
(2) u e X 1 1 oc (R+ +1 ;R m ) is called asymptotically steady if there is a w : R d — > 
R m , so that 

u(t,>)->w mLl c (R d ) 

for almost alH > T. 

Theorem 2. If u £ L°°(R^ +1 ; R m ) is an asymptotically steady and bounded weak 
solution, then w (as in Definitional is a weak solution as well. If u is an EEF 
solution, so is w. 

Proof. Let (f> <E C£°((0, oo) x R d ) be arbitrary. Let supp0 C [0, r]. For any e > 0, 
w(x)<pt (t) + f(, w ( x )) " V^>(i, di 



u(T + t,x)</!>t(i,x) + /(u(T + i,x)) • V0(i,x) dx dt + O^eU-D^oo) 
= 0^11^1100) 

because we can extend </>(• — T, •) e C C ((T, 00) x R d ) smoothly by to a map 
£ C c ((0, 00) x R d ) and use that u is a weak solution. Lcmmanshows that w is 
a weak solution. 

For the proof of the EEF part, replace u,w by rj(u),r](w) and f(u),f(w) by 
ijj(u),ip(w) above. □ 
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